Experimental search for solar axions. 
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A new technique has been used to search for solar axions using a single crystal germanium detector. It exploits 
the coherent conversion of axions into photons when their angle of incidence satisfies a Bragg condition with a 
crystalline plane. The analysis of approximately 1.94 kg.yr of data from the 1 kg DEMOS detector in Sierra 
Grande, Argentina, yields a new laboratory bound on axion-photon coupling of <? aTr < 2.7 x 10 -9 GeV -1 , 
independent of axion mass up to ~ 1 keV. 



1. INTRODUCTION 

The axion is the Nambu-Goldstone boson of the 
broken chiral global U(l) symmetry introduced 
by Peccei and Quinn twenty years ago to dy- 
namically solve the, so called, "strong CP prob- 
lem" of QCD. This prompted many theoretical 
investigations and experimental searches. 

The axion-photon interaction Lagrangian is 
tint = {l/AM)aF av F al3 e nlJa p where a is the 
pseudoscalar axion field, F M „ is the electromag- 
netic field, and 1/M = g a jf is the axion-photon 
coupling. The objective of this experiment is to 
detect solar axions through their coherent Pri- 
makoff conversion (see Fig. [j]) into photons in 
the lattice of a germanium crystal when the in- 
cident angle satisfies the Bragg condition. As 



it turns out [|[, the detection rates in various 
energy windows are correlated with the relative 
orientations of the detector and the sun. This 
correlation results in a temporal structure which 
should be a distinctive, unique signature of the 
axion. We present here the results of a search 
using a 1 kg, ultra-low background germanium 
detector installed in the HIPARSA iron mine in 
Sierra Grande, Argentina at 41° 41' 24" S and 
65° 22' W. A complete description of the exper- 
imental set-up was given earlier by Di Gregorio, 
et al.j| and Abriola, et al. ||. This experiment 
was motivated by earlier papers by Buchmuller 
and Hoogeveen || and by Paschos and Zioutas 
U ; the present technique was originally suggested 
by Zioutas and developed by Creswick, et al. [^J. 
The terrestrial flux of axions from the sun can 



*Supported partially by grants from CONICET, Fun- 
dacion ANTORCHAS and UNSAM 



71-26 




Figure 1. Feynman diagram of the Primakoff con- 
version of axions into photons. 



be approximated by the expression H : 
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10 ) and is dimension- 
1.103 keV, and $ D = 5.95 x 
. The total flux for A = 1 inte- 
grated from to 12 keV is 3.54 x 10 15 cm -2 sec -1 . 
The spectrum is a continuum peaking at about 
4 keV decreasing to a negligible contribution 
above 8 keV. The differential cross section for Pri- 
makoff conversion on an atom with nuclear charge 
Ze is: 
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where q is the momentum transfer, k is the mo- 
mentum of the incoming axion, and r is the 
screening length of the atom in the lattice. For 
germanium, a — Z 2 ah 2 c 2 g ail /&n = 1.15 x 
10~ 44 cm 2 when g ail = 10~ 8 GeV -1 , or equiva- 
lently A = 1. 

For light axions the Primakoff process in a peri- 
odic lattice is coherent when the Bragg condition 
(2d sin 9 = nX) is satisfied, that is when q trans- 
ferred to the crystal is a reciprocal lattice vector 
G = 2n(h, k, l)/a . Here a a is the size of the con- 
ventional cubic cell, and h, k, and / are integers. 



It was shown that the rate of conversion of ax- 
ions with energy E when the sun is in the direc- 
tion k, N(k,E), can be written [||: 



N(k,E) = 2hc 
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where V is the volume of the crystal, v c is the 
volume of a unit cell, S(G) is the structure func- 
tion for germanium, and d$/dE is evaluated at 
the axion energy of h~c\G\ 2 /2k ■ G. The structure 
function for germanium is: 
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Note that in (|^) the coherent conversion of axions 
occurs only for a particular axion energy given the 
position of the sun, fc, and reciprocal lattice vec- 
tor G. However, the detector has a finite energy 
resolution; for the detector in Sierra Grande it is 
1 keV FWHM at 10 keV. We take this into ac- 
count by smoothing N(k, E) with a Gaussian of 
the appropriate width. Finally, we take the rel- 
evant part of the energy spectrum, in this case 
from the threshold energy of 4 keV up to 8 keV 
(which is just below the X-rays at 10 keV), and 
calculate the total rate of conversion in windows 
of width AE, typically 0.5 keV, 
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where E a = hc\G\ /(2k ■ G) and erf(x) = 
Jq e _t dt is the error function. In equation 
(||) we have neglected the angular size of the core 
of the sun and the mass of the axion which is jus- 
tified when m a c 2 is small compared to the core 
temperature of the sun ||, i.e., up to a few keV. 

The theoretical axion detection rate for this de- 
tector, calculated with equation (|J), is shown in 



Figure g. The position of the sun is computed at 
any instant in time using the U.S. Naval Observa- 
tory Subroutines (NOVAS) ||. The pronounced 
variation in R(k,E) as a function of time invites 
the data to be analyzed with the correlation func- 
tion: 
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Figure 2. A typical axion-photon conversion rate, 
R{t,E), for various energy bands. The experi- 
mental energy resolution FWHM=1.0 keV at 10 
ke V was used. 



X = ]T [R(k,E) - (R(E))] n(U) , (6) 

i=l 

where R(ti, E) is the smooth shape of the theoret- 
ical rate at the instant of time, ti, (R(E)) is the 
average rate over a finite time interval, and n(ti) 
is the number of events at ti in a time interval 
At, usually or 1. The choice for the weighting- 
function W(t,E) = R[k(t),E] - (R{E)) is moti- 
vated by the requirement that any constant back- 
ground average to zero in x, whereas a counting 
rate which follows R[k{t),E] increases x- 

The number of counts at time, t, in the interval 
At is assumed to be due in part to axions and in 
part to background governed by a Poisson process 
with mean: 

(n(t)) = [XR(t,E) + b(E)} At, (7) 



where b(E) is constant in time. 
The average value of x 1S then, 

( X ) = Y^muE) - (R(E))][XR(ti,E) + 

i 

b(E)} At 
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i 

We can add and subtract the constant quantity 
X(R(E)) to the second factor in eq. (||). Any time 
independent contributions multiplied by W(t, E) 
in cq. (^|), and summed over time, will vanish. 
Accordingly, taking the limit as At — > 0, we ob- 
tain: 

( X (A)) = X f W 2 (t,E)dt. (9) 
Jo 

The expected uncertainty in x, (Ax) 2 = (x 2 ) — 
(x) 2 , is given by, 

(A X ) 2 = ^^W^mWfaE) x 

i 3 
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i 

(10) 

where the square bracket is (An(ti)) 2 , which in 
Poisson statistics is (n(ti)). Accordingly, 

(Ax) 2 - ^W^mniU)). (11) 

i 

By (Q) we have: 

(A X ) 2 = J2 w2 ^ E ) [*R(.U,E) + b(E)} At 

i 

= ^2 W 2 (U, E) {X [R(U, E) - (R(E))} 

i 

+ X(R(E)) + b(E)} At, (12) 
which in the limit At — > becomes, 

(Ax) 2 = A f W 3 {t,E)dt + 
Jo 

R T {E) [ W 2 {t,E)dt. (13) 
Jo 



The quantity Rr(t, E) is the average total count- 
ing rate, including both axion conversions and 
background. 

The data are separately analyzed in energy 
bins, AEk, fixed by the detector resolution 
(FWHM - 1 keV in this case). The likelihood 
function is then constructed: 



L(A) = J] exp 



~(Xk ~ {Xk)f 
2(A Xfc ) 2 



(14) 



to place a bound on the axion interaction rate, 
the data must be analyzed for many azimuthal 
orientations of the crystal, and the weakest bound 
selected. The results of these calculations for 
707 days of data in the energy range from 4 to 
8 keV in 0.5 keV intervals are shown in Figure ||. 



To an excellent approximation (Ax^) 2 is domi- 
nated by background. Maximizing the likelihood 
function, the most probable value of A is given by 



An 



where, 
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and the width of the likelihood function is given 
by, 



-1/2 



(17) 



We note that A k is proportional to the time of 
the experiment, so that o\ decreases as T 1 / 2 . 
The background scales with the detector mass, 
while Ak scales as the square of the detector mass, 
therefore a\ decreases as (M^T) -1 / 2 . 

We have carried out extensive Monte-Carlo test 
of this method of analysis. As a test of our anal- 
ysis, typical results for the likelihood function for 
the cases A = (no axions) and A = 0.003 were 
calculated with realistic backgrounds for a detec- 
tor operating with the same mass, energy reso- 
lution, and threshold as the DEMOS detector at 
the latitude and longitude of Sierra Grande for 
one year. It is clear from this calculation that 
the correlation function analysis is consistent and 
quite sensitive to the presence of a variation in the 
counting range due to solar axions with a signal 
to noise ratio less than 1%. 

The vertical axis of our detector is the (100) 
crystalline axis. The orientation of the (010) and 
(001) axes are unknown at this time. Therefore, 
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Figure 3. Values of A calculated from the 707 
days of data as a function of the azimuthal angle 
6. The error bars are la. 



The most conservative upper bound on A, or 
equivalently gayy, is found by taking for <f> the 
angle at which A is maximum. This yields an up- 
per bound on the axion-photon coupling constant 
g ail < 2.7 x 10~ 9 GeV" 1 at the 95% confidence 
level. 

In Figure ^ we show the area of the axion mass 
- coupling constant plane excluded by this result 
along with results of earlier work. 

While this bound is interesting because it is a 
laboratory constraint, it does not challenge the 
bound placed by Raffelt §, g ail < lO" 10 GeV" 1 
based on the helium burning rate in low mass 
stars. A coupling constant g ail ~ 10 -9 GeV -1 
would imply axion emission rates 100 times higher 
than the stellar bound, and a significantly differ- 
ent concept of stellar evolution. 
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Figure 4. Exclusion plots on the g Q77 us. axion- 
mass plane. The curve to the left are from 
Ref. /[7^/. TAe letters indicate the following he- 
lium gas pressures in the conversion region of a 
strong magnetic field: (a) vacuum, (b) 55 Torr, 
and (c) 100 Torr. 



This experiment can be considerably improved 
by using a large number of smaller p-type ger- 
manium detectors with known orientations of 
the (010) axes, with energy thresholds below 
2 keV, and energy resolutions corresponding to 
FWHM « 0.5 keV. This is being proposed at 
this time. Another collaboration could also op- 
erate the COSME experiment in the new Univer- 
sity of Zaragoza underground laboratory in the 
Canfranc tunnel at 42° 48' N and 0° 31' W. The 
COSME detector is a 0.25 kg crystal having an 
energy threshold of ~1.8 keV and a resolution 
of 0.5 keV FWHM at 10 keV. A positive result in 
this northern hemisphere experiment over a wider 
energy range should have a very different tempo- 
ral pattern from that of the Sierra Grande exper- 
iment, but should yield the same value of A. 

One of the USC/PNNL twin detectors is cur- 
rently operating in the Baksan Neutrino Obser- 
vatory in Russia at 660 mwe, and if moved to 



a location with greater overburden could also be 
used to acquire meaningful data on solar axion 
rates. 

The analysis presented here allows one to legiti- 
mately combine the results of a number of experi- 
ments. Accordingly, the results from a large num- 
ber of experiments located throughout the world 
can be combined to yield results equivalent to a 
single large experiment. 
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